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Abstract 

This  paper  explores  further  the  use  of  Interelectron  coordinates 
in  constructing  atomic  wave  functions.   A  simple  method  is  developed 
for  constructing  wave  functions  of  this  type  which  yields  surprisingly 
good  values  for  the  energy  of  a  variety  of  atomic  systems,  in  zero  order. 
The  hamiltonian  for  the  system  is  split  into  an  unperturbed  part  which  is 
separable  and  which  contains  the  interelectron  potentials  as  well  as  the 
electron-nucleus  potentials  and  into  a  perturbing  term  which  is  always 
finite  and  which  vanishes  whenever  an  electron  is  far  from  the  nucleus. 
These  zero-order  energies  corresponding  to  this  splitting  of  the  hamil- 
tonian are  at  least  an  order  of  magnitude  better  for  the  light  atoms 
than  the  energies  given  by  the  usual  Thomas-Fermi  theory  and  are  consi- 
derably better  than  the  energies  calculated  with  hydrogenic  functions 
alone  in  first-order. 
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1.   Introduction 

Several  years  ago,  Pluvinage  '-  -'  illustrated,  for  helium,  a  method 
in  which  interelectron  coordinates  could  be  used  to  separate  partially  the 
Schrodinger  vave  equation  so  that  both  the  interelectron  potentials  and  the 
electron-nucleus  potentials  were  taken  into  the  unperturbed  part  of  the 
hamiltoniaa.  In  the  first  paper  of  this  series'--"  Pluvinage's  approach 
was  modified  and  an  approximation  technique  was  developed  so  that  varia- 
tional values  of  the  energy  of  systems  beyond  helium  could  be  accurately 
calculated.   The  groimd  state  energy  of  lithium  was  calculated  and  the 
result  obtained  indicated  that  Pluvinage's  approach  held  excellent  promise 
for  other  systems. 

In  this  paper,  the  method  of  Pluvinage  is  modified  in  a  different  way 
in  order  to  obtain  relatively  simple  zero-order  wave  functions  which  include 
correlation.   The  unperturbed  hamiltonian  is  separable  and  the  resulting 
differential  equations  are  hydrogenic  in  character.   Thus,  the  ground 
state  or  the  excited  state  energy  of  a  variety  of  atomic  systems  can  be 
explicitly  written  down.   The  resiilts  obtained  are  considerably  better  than 
those  obtained  using  hydrogenic  functions  in  zero-  or  in  first-order,  or 
using  the  Thomas-Fermi  method.   The  zero-order  energies  we  calculate  are 
surprisingly  close  to  variational  energies  calculated  using  the  wave  func- 
tions without  explicit  r.  .  dependence. 

The  motivation  behind  this  paper  was  to  obtain  a  separation  of  the 
hamiltonian  for  atomic  systems  which  leads  to  a  solvable  \inperturbed  pro- 
blem and  to  a  perturbing  term  which  is  truly  small.   In  addition,  of  the 
many  possible  solutions  to  the  unperturbed  hamiltonian  only  those  are  chosen 
which  display  the  physical  behavior  expected  of  the  system. 


The  material  in  [2]  on  the  two  electron  system  had  previously  been  done 
but  not  published  by  M.  Kelly  and  L.  Spruch.   Reference  to  their  contribu- 
tion was  inadvertently  omitted  in  our  first  paper. 


The  actual  method  adopted  was  as  follows .   A  full  set  of  dependent 
Interelectron  coordinates  is  introduced  in  addition  to  the  usual  independent 
coordinates  of  the  electron  with  respect  to  the  nucleus.   We  shall  refer  to 
these  electron  coordinates  with  respect  to  the  nucleus  simply  as  nuclear  coor- 
dinates from  now  on.   With  these  variables,  the  hamlltonian  can  be  split  into 
an  unperturbed  part  and  a  perturbing  part.   The  unperturbed  hamiltonian  is  it- 
self separable  in  the  various  nuclear  and  interelectron  coordinates.   It  con- 
tains all  of  the  interactions  present  in  the  system  and  is  exactly  solvable. 
The  perturbing  hamlltonian  differs  from  that  used  previously"-  -1 '  L  J  in  that  it 
vanishes  whenever  an  electron  is  far  from  the  nucleus  while  still  remaining 
finite  when  any  coordinate  vanishes.   This  perturbing  term  is  small  compared  to 
the  Coulomb  interactions  in  the  neighborhood  of  their  singularities  and  is  com- 
parable to  the  Coulomb  potentials  at  large  distances  from  the  atom  where  the 
Coulomb  potentials  and  the  perturbing  term  are  both  small.  We  then  hope  that 
the  perturbing  term,  chosen  this  way,  would  be  small  throughout  all  of  space. 
The  zero-order  wave  function  is  a  product  of  the  individual  nuclear  and  inter- 
electron wave  functions  while  the  zero-order  energy  is  a  sum  of  the  indivi- 
dual nuclear  and  interelectron  energies.   To  specify  in  detail  the  zero-order 
solutions,  suitable  boundary  conditions  must  be  introduced.   As  usual,  the 
part  of  the  wave  function  expresses  in  the  nuclear  coordinates  must  vanish 
ar  infinity  and  be  finite  at  the  nucleus.   This  leads  to  the  bound  state 
hydrogenic  functions  for  this  part.   We  impose  boimdary  conditions  on  the  inter- 
electron wave  functions  to  account  for  the  physical  behavior  expected  of  the 
system;  namely,  that  when  an  individual  electron  is  far  from  the  nucleus,  the 
nuclear  charge  is  shielded  by  the  remaining  electrons.   The  part  of  the  wave 
function  containing  the  interelectron  coordinates  is  required  to  account  for 


this  shielding  at  large  nuclear  distances.  This  iiniquely  determines  this 
part  of  the  wave  function  and  leads  immediately  to  appropriate  zero-ordar 
energies . 

In  Sec.  2  of  this  paper,  the  general  form  of  the  hamiltonian  is  pre- 
sented and  discussed  vhen  the  interelectron  coordinates  are  introduced  in 
addition  to  the  radial  and  angular  coordinates  generally  used.   In  Sec.  3/ 
the  separation  of  the  hamiltonian  into  a  perturbed  and  an  unperturbed  part 
is  presented.   The  boundary  conditions  to  be  used  on  the  radial  and  inter- 
electron parts  of  the  vave  function  are  given  in  Sees,  k   and  5-   The  solu- 
tion of  the  unperturbed  vave  equation  londer  the  boimdary  conditions  imposed 
is  easily  carried  out  and  the  zero-order  energy  is  obtained  immediately. 
In  Sec.  6,    application  is  made  to  various  charged,  neutral  and  excited  atoms 
to  illustrate  the  method  and  its  results.   The  results  are  discussed  in 
Sec.  7  and  comparison  is  made  there  with  other  methods  for  calculating  atomic 
energies. 

2.   The  Hamiltonian  with  Interelectron  Coordinates  , 

The  hamiltonian  for  the  Schrodinger  equation  is  written  in  atomic  units 
as 

1  11       1^,1   ij 

Z  is  the  nuclear  charge  and  the  summations  run  over  the  N  electrons 
present  within  the  atom.   The  interelectron  potentials,  l/r.  .,  prevent  the 
exact  solution  of  this  Schrodinger  equation.   The  fact  that  these  potentials 
are  singular  and  numerous  prevents  any  great  success  in  treating  them  as 
perturbations.   This  difficulty  can  be  overcome  to  a  great  extent  by  intro- 


ducing  the  interelectron  coordinates  explicitly. 

The  usual  coordinates  taken  to  represent  the  hamiltonian  (2.1)  for  an 
atomic  system  are  the  JN  radial  and  angular  coordinates  r.^G.^(()..  These  are 
the  nuclear  coordinates.  If  we  nov  introduce  additional  interelectron  coor- 
dinates then  ve  have  more  coordinates  to  describe  the  system  than  degrees  of 
freedom  present  within  the  system.  This  difficulty  is  overcome  by  consider- 
ing the  r.  .  as  convenient  dependent  variables  which  are  expressible  in  terms 
of  the  true  independent  variables  chosen  for  the  system.   Thus 


r.  +  r.  +  2r.r.(cos  0.  cos  0.  +  sin  0.  sin  0.  cos(d)  _(t)  )) 


=    V(x,-  X.)  +  (y.-  y.  )  +  (z.-  z.  )   . 
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J   1 


(2.2) 


In  order  to  express   the   hamiltonian   of  the   system  in   terms   of  the   full   set 
of  variables^    dependent  plus   independent^    the   usual  rules    for  differentiating 
are  used.      The   result   is 


H 
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where  0.  and  (p.  are  unit  vectors  given  by 

1     ^1 


G.  =  cos  9.  cos  6.  i  +  cos  9  sin  6.  j  -  sin  9. 
^.  =  -  sin  (b.  i  +  cos  (f).  j. 


(2.4) 


Note  that  the  introduction  of  the  dependent  r. .  does  not  destroy  the 
orthogonality  among  the  r.,9.  and  (j).  and,  consequently,  there  are  no  terms 
coupling  these  variables  directly.   The  lack  of  orthogonality  between  these 
variables  and  the  r. .,  and  among  the  r. .,  account  for  the  coupling  terms  which 
are  present. 

Extending  Pluvinage's  idea,  we  can  see  from  (2.3)  how  the  total  potential 

energy  coiold  be  included  in  zero  order.   The  first  two  summations  on  the  left 

might  be  taken  to  represent  the  unperturbed  hamiltonian  (we  shall  not  do 

exactly  this).   The  corresponding  zero-order  Schrodinger  equation  is  then 

separable  in  the  r.,9.,(t).  and  r.  .  coordinates.   The  zero-order  solution  is  a 

product  of  the  individual  solutions  for  each  of  the  coordinates,  dependent  and 

independent,  while  the  zero-order  energy  is  a  sum  of  the  zero-order  energies 

obtained  for  each  coordinate  when  the  proper  boundary  conditions  are  imposed. 

The  fact  that  the  r. .  are  dependent  variables  makes  it  difficult  to  carry  out 

integrations  involving  this  solution  over  all  space  to  determine,  say,  the 

normalizing  constant.   The  element  of  volume  is  IT  r.  sin  6.  dr.  d9.  dd).  . 
•^  .   1      1   1   1   ^1 

1 
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The  r.  .  -vd-thin  the  zero-order  vave  function  must  first  be  expressed  in  terms 
of  the  radial  and  angular  variables  by  Eq.  (2.2)  and  then  the  integrations  per- 
formed.  For  any  complicated  system_,  such  integrations  must  invariably  be  done 
approximately  and  one  particular  method  of  doing  so  accurately  was  described 
in  1.   This  emphasizes  the   importance  of  obtaining  the  best  zero-order  solu- 
tion possible  in  the  hope  that  accurate  information  can  be  obtained  by  examin- 
ing just  the  zero-order  case. 

3.    Separation  of  the  Hpimiltonian 

A  possible  method  of  separating  the  hamiltonian,  once  the  interelectron 

coordinates  are  introduced^  was  suggested  in  the  previous  section.   For  reasons 

given  below^  the  actual  separation  chosen  is  somewhat  different.   Instead^  we 

set  H  =  H  +  H'  where 
o 


H  = 
o 


1 
2 
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The  choice  given  above  differs  from  that  suggested  in  the  previous 

section  by  the  addition  of  second  derivative  terms,  -  -^r-  b.  .  S  /Sr.  .  ,  in  the 

2  ij   '   ij 

unperturbed  hamiltonian  and  their  subsequent  subtraction  in  the  perturbing 
part.  We  use  the  freedom  which  we  have  gained  by  introducing  the  arbitrary 
constants  b . .  to  make  the  perturbation  vanish  in  the  asymptotic  regions  of 
space. 

The  zero-order  wave  function  is  then 


%-    TT  ($i(-i^e.,^.)   IT   u  (r  )}.  (3.3) 

i^  j>i'^ 

The  ^.    are   the   hydrogenic  wave   functions'--'    for   an   atom  of  charge   Z 

$i(r.,9.,(|).)    =  R^   ^    (r.)Y,  ^  (9.,(|).  )  (3.4) 

i   i  i    i 

where  R  represents  the  radial  ei gen functions  while  Y  represents  the  angular 

spherical  harmonics,   n, ^  and  m  represent  the  appropriate  orbital  quantum 
numbers. 
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The  u. .  satisfy  the  following  equation 


d  u.  .  „  du.  .   u.  . 

(1  +  b.  .)  |i  -  ^-  -^  +  .J^  =  E.  .u.  .  .                    (5.5) 

ij'    2  r. .  dr. .   r. .    ij  ij                     ^^  ^' 

dr^^  ij  ij    ij 


The  solution  of  this  equation  depends  on  vhat  we  choose  for  the  asymp- 
totic behavior  of  the  u.  ..   We  shall  insist  that  \|f  be  square  integrable.   In 

ij  o     ^        ° 

section  5,  we  shall  see  that  in  order  to  accoixnt  for  the  shielding  of  the  i 
electron,  we  require  an  exponentially  increasing  solution  of  (5-5)^  i.e., 
u.  . -^  exp(q_^  .r.  .),  and  that  the  shielding  requirement  enables  us  to  determine 
the  q. ..   The  solution  of  (5«5)  is  not  an  eigenvalue  problem.   However,  in 
terms  of  the  still  undetermined  b. .  and  q. .,  the  energy  is  given  by  (see 


Appendix  I ) 


E.  .  =  -  (1  +  b.  .)  q.^  (5.6) 


The  scheme  for  determining  the  b. .  will  be  discussed  in  section  h.      Thus 
we  shall  be  able  to  determine  the  unperturbed  energy  knowing  nothing  but  the 
asymptotic  form  of  the  u. .. 

The  proper  antisymmetrization  of  the  zero-order  wave  function  in  Eq. 
(5.5)  is  a  simple  generalization  of  well-known  procedures.   As  expected,  the 
antisymmetrization  requires  that  no  two  electrons  have  all  four  quantum  num- 
bers (spin  included)  the  same.   At  most  two  electrons,  each  of  opposite  spin, 
can  occupy  the  same  atomic  orbit.   The  only  way  in  which  we  will  make  use  of 
the  antisymmetrization  is  by  applying  the  Paul!  principle  in  filling  these 
atomic  orbits.   Otherwise  the  zero-order  wave  function  will  be  used  in  its 
unsymmetrized  form.   One  result  Ol  this  is  that  the  unperturbed  hamlltonian 


vill  not  be  symmetrical  in  all  of  the  electrons.   The  antisymmetrization  is 
discussed  briefly  in  Appendix  II, 


h.        Uetermlnation  of  the  b. . 
ij 

The  total  unperturbed  energy  of  the  system  depends  on  the  choice  of 
two  arbitrary  parameters,  b. .  and  q. ..   The  latter  are  fixed,  as  has  been 
indicated,  by  the  asymptotic  behavior  of  the  vave  function.   The  former  con- 
stants are  at  our  disposal  and  we  choose  them  to  minimize  the  effects  of  H' . 
The  potential  terms  in  our  hamlltonian 


i    i       j  ?t  k    jk 


are  singular  at  the  zeros  of  r.  and  r. ..   Since  the  potential  energy  is  com- 
parable to  the  kinetic  energy,  V  cannot  be  thought  of  as  a  small  perturbation. 
The  introduction  of  the  interelectron  coordinates  allows  both  the  nuclear  and 
the  interelectron  potentials  to  be  taken  into  the  unperturbed  hamiltonian,  H  . 
The  H  given  by  Eq.  (5-1)  contains  all  of  these  potentials  and  allows  their 
behavior  at  the  singularities  to  be  treated  accurately. 

Since  all  the  Coulomb  Interactions  have  now  been  included  in  the  uhper- 
turbed  hamlltonian,  we  can  hope  to  obtain  accurate  zero-order  results  if  the 
perturbing  hamlltonian  is  small  compared  to  V.  At  the  zeros  of  r.,r  ,  we 

1   JK 

note   from  Eqs.    (5.2)    and   (5.5)   that  H'lir     is   always  finite.      (This   is   obvious 
for  the   zero  angular  momentum  cases,   where  the   angijlar  derivatives   are   zero. 
i?'or  other  angular  momenta  the  behavior  of  the  R(r   )    cancels   the  r.      in  the 
angular  derivative   terms   near  zero).      As   already  mentioned,    V>!'     becomes 
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infinite  at  the  points  in  question  and  the  effects  of  the  perturbing  hamil- 
tonian  are  thus  small  compared  to  the  zero-order  potentials  in  the  regions 
of  space  near  the  nucleus  and  near  each  electron^  independent  of  the  choice 
of  any  arbitrary  constants. 

The  problem  is  more  complex,  however,  in  vhat  ve  will  refer  to  as  the 
asymptotic  region  of  space,  where  an  electron  is  far  from  the  nucleus  and 
its  companion  electrons.   In  this  region  we  will,  of  course,  want  to  compare 
H'  with  V.   Then,  by  adjusting  the  b. .,  we  will  make  H'  vanish  in  a  suitable 
manner.   The  first  thing  to  be  done  is  to  decide  what  region  of  space  shall 
comprise  the  asymptotic  region.   After  this,  we  will  be  able  to  detail  the 
suitable  manner  in  which  H'  must  vanish. 

In  order  to  determine  the  asymptotic  region  of  space  for  each  electron 
we  picture  the  atom  as  being  constructed  by  bringing  the  electrons  in  from 
infinity,  filling  the  atomic  shells  (determined  by  the  principal  quantum 
number,  n)   consecutively,  starting  from  the  innermost.  In  eux   approximation, 
an  electron  in  an  inner  shell  is  not  affected  by  any  outer  shell,  but  it  is 
by  the  electrons  in  lower  shells  or  in  the  same  shell.   Such  a  picture  is 
simple  to  apply  and  allows  a  unique  determination  of  the  shielding  require- 
ment imposed  as  a  boundary  condition  upon  the  interelectron  wave  functions 
in  the  next  section.   Its  applicability  is  discussed  in  the  last  section  of 
this  paper. 

Consider  an  electron,  the  i   ,  which  occupies  the  n  shell.   The  atom 
Is  built  up  to  and  including  this  shell.   The  asymptotic  region  for  this 
electron  is  created  by  moving  the  electron  outward  from  the  vicinity  of  the 
incomplete  atom.   The  Coulomb  potential  acting  upon  the  electron  becomes 
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V.  =  -  ^^  + 
1     r.     .■'— ; — .   r. 

1    ^  =^   J        ij 


Z    ^ 


-Z  +  P.-  1 

1 


(^^.2) 


where  P.  is  the  total  number  of  electrons  in  all  shells  up  to  and  including 
the  n. .   in  the  perturbing  hamiltonian^  the  portion  containing  the  coordi- 
nate of  the  i   particle  is 


1    [— 


H 


J  ?t  1 


r.  T.  .     -.c 

_J^ ij^   o. 
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In  the  asymptotic  region  the  individual  vectors  r.  and  r. .  line  up.   The  dis- 
tance derivatives  have  the  effect  of  multiplying  t  by  the  constants  appeax- 
ing  in  the  exponents  of  R  and  u,  vhile  the  angular  derivative  terms  become 
negligible  because  of  the  extra  factor  of  r.  in  the  denominator.   Then  the 
operator  H.'   has  the  following  approximate  numerical  value 


e: 


-b. 


ij 


%.i  -  '^.l 


n.   2   .4-7- 


t/j,i 


^ij  'hi 


+  a(r:^)  {k.k) 


vhere  ve  have  again  assumed  that  u.  .  — >  exp(q_^  .r.  .j  in  the  asymptotic  region. 
The  sum  is  over  all  shells  up  to  and  including  the  n.  shell.   This  is  indica- 


ted by  the  prime  on  the  summation  sign. 
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The  perturbation  caused  by  the  i   electron  thus  behaves  like  a  con- 
stant in  the  asymptotic  region,  while  the  corresponding  Coulomb  potential 
found  in  the  imperturbed  term  vanishes  like  l/r. .   This  is  not  too  disturbing 
since  HI  ^     will  vanish  because  of  the  exponentially  decreasing  behavior  of 
i|f  .   Certainly,  however,  an  improvement  is  expected  if  the  b.  .  we  introduced 
are  adjusted  so  that  H!  itself  will  vanish  at  least  as  strongly  as  V.  in  the 
asymptotic  region.   Therefore  we  take 


/     -i3. ,  q.  .  -  q.  .  —  +7^      /  q.  .  q.. 


0  *   1 


ij  ^0    ^xj  "i   -      I  ^  j 


=0.  (^.5) 


This  is  the  prescription  we  use  to  insure  that  the  perturbation  term 
is  small  in  the  asymptotic  region,   it  is,  of  course,  applied  to  all  the 
various  shells  comprising  the  atom.   Note  that  only  n.  enters  into  the  formula, 
so  that  all  electrons  in  the  same  shell,  but  with  different  m, I  values,  are 
treated  identically.   The  q's  are  determined  by  the  shielding  discussed  in  the 
next  section. 

We  now  have  completed  the  description  of  the  separation  of  the  hamlltonian 
and  the  determination  of  the  b. ..   The  perturbation  introduced  is  small  com- 
pared  to  the  Coulomb  potentials  near  the  origins  of  all  the  coordinates.   In 
the  asymptotic  region  about  each  electron,  the  effect  of  the  perturbation, 
although  comparable  to  the  Coulomb  potential,  decreases  as  l/r.   It  is  then 
reasonable  to  expect  that  the  effect  of  the  perturbation  will  be  small  through- 
out nil  of  the  physical  space  of  each  electron,  and  that  the  zero-order  wave 
function,  Eq.  (5.5),  should  accurately  represent  the  system. 
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The  asymptotic  behavior  of  H!  for  the  interelectron  wave  function  used. 
in  I  is  worthy  of  note.   The  q's  were  all  taken  as  zero  in  that  paper  and  all 
the  interelectron  wave  functions  became  the  same  Bessel-Clifford  function. 
The  corresponding  HI  thus  vanish  asymptotically.  However^  investigation  of 
the  Bessel-Clifford  function  shows  that 


K   *BC  ~  ^(-I'^')  ^BC 


In  addition^  the  shielding  condition  we  impose  in  the  next  section  cannot  be 
met  by  this  wave  function.  Nevertheless,  the  approach  used  in  I  resembles 
in  part  the  present  method  and  the  success  found  with  the  Bessel-Clifford 
function  partially  validates  the  method  used  here. 

5-   The  Shielding  Requirement 

The  zero-order  wave  function  has  been  given  in  Eq.  (5*5) •   In  order  to 
determine  the  wave  functions  J  and  u  completely,  suitable  boundary  conditions 
must  be  imposed.   Actually  in  writing  J  as  given  in  Eq.  (5.^)  we  have  impli- 
citly assumed  one  set  of  boimdary  conditions.   This  is  that  the  radial  wave 
functions  are  finite  at  the  origin  and  vanish  at  infinity.   In  fact,  they 
must  vanish  so  strongly  at  infinity  that  the  total  wave  function  t  is  quad- 
ratically  integrable.  Despite  the  exponentially  Increasing  behavior  which 

we  will  demand  of  the  u.  .,  the  solution  for  \lf  satisfies  this  boundary  condl- 

ij  o 

tion. 

What  boundary  condition  shall  we  impose  upon  the  u.  .?  The  choice,  we 
have  indicated  earlier,  must  express  the  shielding  of  the  nuclear  charge  which 
occurs  as  one  electron  is  moved  away  from  the  nucleus. 
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Let  us  examine  the  solution  of  the  Schrodlnger  equation  vhen  the  i 
electron  in  the  n.  shell  is  in  its  asymptotic  region  as  described  in  the 
previous  section.  Mathematically  we  do  this  by  requiring  that  r  /r.  be  very 
small,  where  the  subscript  h  designates  any  electron  in  the  same  shell  or  in 
an  inner  shell  compared  to  i,  while  requiring  that  r  /r.  be  very  large  where 
the  subscript  k  represents  any  electron  in  the  outer  shells.   The  cloud  of 

electrons  located  in  shells  outside  of  the  n.  shell  will  not  affect  the 

1 

potential  of  the  i   electron,  assuming  that  these  electrons  are  uniformly 
distributed.  We  specifically  assume  that  such  is  the  case,  although  there 
are  errors  involved  in  this  assumption  due  to  the  overlap  of  the  wave  func- 
tions corresponding  to  different  shells.   The  inner  electrons,  h,  however, 
will  shield  the  nuclear  charge,  changing  its  effective  value  to  Z  -  P.+  1. 

That  part  of  the  Schrodinger  equation  corresponding  to  the  i   electron 
now  separates  out.   If  we  impose  the  usual  boundary  conditions  on  the  solution 
of  this  part  of  the  Schrodinger  equation,  that  it  be  finite  at  small  r  and 
vanish  at  infinity,  the  corresponding  wave  function  is  a  hydrogenic  wave 
f\mction  bound  in  some  orbit  to  a  charge,  Z  -  P.+  1.   Under  the  conditions 

that  r,  /r  must  be  small  while  r,  /r.  is  large,  the  i   electron  does  not  have 
h'  1  k'  1 

any  possibility  of  interacting  with  the  other  electrons,  and  it  must  remain 
in  its  original  orbit  in  the  n.  shell.   Therefore,  in  the  asymptotic  region 
of  the  i   electron,  the  zero-order  wave  function  must  have  a  separable  part 
which  represents  the  i   electron  as  moving  in  its  original  shell  in  the  field 
of  a  nuclear  charge  Z  -  P.+l.   This  requirement  cannot  be  met  in  every  detail. 
However,  it  is  easily  satisfied  as  far  as  the  exponential  behavior  of  *  is 
concerned  by  imposing  the  following  boundary  condition  on  u.  .  at  large  dls- 
tances: 
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P. 
$i(r.)   _  J\_      ^i/^ij^  ^    exp[-(Z-P.+l)/n.]r.  ,  (5.I) 


Thus 

P. 
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p.- 1 
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•  ^  .   ij     n. 


since  J.(r.j  ^-'  exp[-(z/n.  jr.]  as  r.  — >  00  . 

In  view  of  the  equivalent  shielding  by  all  the  inner  electrons,  we  take 


all  of  the  q. .  as  equal.   That  is 


1 

1-  ■  = 

) 

n  .   <  n. 

IJ 

n. 

1 

J  -     1 

(5.2) 


Since  the  outer  electrons  do  not  shield  the  i   ,  we  need  not  impose  any 
requirement  on  the  q.  .,  n.  >  n.  .   The  arg\:ment  given  above  is,  of  course, 
not  a  proof.  However,  it  does  establish  the  reasonableness  of  our  approach. 
The  reason,  for  instance,  why  we  do  not  define  the  asymptotic  region  for  the 
i   electron  by  moving  this  electron  out  past  all  of  the  other  electrons 
becomes  obvious  from  our  argument.   In  passing  the  i   electron  through  the 
cloud  of  outer  electrons,  the  validity  of  the  original  orbital  assignment 
will  be  destroyed  because  of  the  electronic  interactions.  We  cannot  say  what 
precise  orbit  the  i   electron  will  now  occupy  around  the  shielded  nucleus 
of  charge  Z  -  W  +  1.   Again,  the  similar  exponential  behavior  of  different 
orbital  wave  functions  representing  the  same  shell  explains  why  the  treat- 
ment is  identical  for  electrons  in  different  orbits  within  the  same  shell. 
Note  that  only  the  principal  quantum  nimiber  appears  in  our  equations. 
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Once  we  have  determined  the  value  of  the  q's,  the  equation  in  the  last 
section  which  expresses  the  condition  that  H!  vanish  asymptotically  can  be 
solved.   Following  the  above  arguments  all  the  b. .  are  taken  as  equal  when 

-'-  (J 

n.  is  less  than  or  equal  to  n..   Then  Eq.  (^.5)  yields 

p. 

-(1  +  b.^)  =  (z  -  ^)  ,    '^j  ^^i  •  (5.3) 


By  starting  with  the  innermost  shell  and  progressively  filling  the  outer 
shells,  a  set  of  numbers  q. .  and  b. .  are  easily  generated.   These  q's  and  b's 
complete  the  description  of  H  and  \|f  .   Useful  properties  of  the  system  can 
now  be  obtained  by  investigating  H  or  \|f  directly.   The  zero-order  energy, 
for  example,  is  just  the  siun  of  all  the  individual  nuclear  and  interelectron 
energies.   From  Eqs.  (3«6),  (5.2),  (5«3)  and  the  properties  of  the  ^'s,  we 
have 


E=/    E.  +  /   /    E.. 
o    ^—       1    ^>V   J^ 


=  -|  H  4-  ^H  (1  +  ^j^  %j  •  (5.^) 

in      i  >  j         J-J    -LJ 

1         ^ 

Various  applications  of  this  equation  will  now  be  discussed. 

6.   Applications 

a.   Two-electron  atoms 

Helium  and  heliimi-like  ions  consist  of  two  electrons  both  in  the  lowest 
orbit  surrounding  a  nucleus  of  charge  Z.   These  are  the  simplest  many-electron 
systems  which  are  encountered.   Only  the  ground  state  energy  is  considered  now. 
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The  zero-order  wave  function  is  a  product  of  two  ground  state  hydro- 
genic  wave  functions  times  the  interelectron  wave  function.   That  is 

-Zr.   -Zr 

to  =  e   ^  e   ^  u(r^2J.  (6.1) 

When  one  electron  is  far  from  the  nucleus  it  moves  in  the  ground  state  of  a 
nucleus  of  shielded  charge^  Z-1.   Its  nuclear  coordinate^  say,  r^,  becomes 
veiy  close  to  the  value  of  r,  p.   The  exponential  from  of  u  at  large  values 
will  be  exp[qr  J  and  in  order  that  t  display  the  proper  shielding  at  large 
r,  we  demand 

-Zr^  -Zr^   qr^2      "^^1  -^^-1^^2  ,,,, 

e     e     e       --^  e     e  <  [o.2) 

Thus 

q  =  1  (6.3) 

in  agreement  with  Eq.  (5.2). 

The  perturbing  hami Itoni an,  H'(r  ,rp)  for  the  second  electron  is  given 
by  Eq.  (J+.3) 

m  t  =  (  b  -^    -    "^£121    ^ —  )  ^  .  (6.4) 

2  ^o    ^  ^^J     r^  r^2  ^^^^^}      ° 

->      — ♦• 

At  large  v  ,    the  vectors  r  and  r   line  up  and  H'  approaches  the  numerical 

value 

H^  =  b  q^  +  Z  q  (6.5) 

In  order  to  have  H'  vanish  asymptotically,  we  require 

b  =  -  Z  ■  {^'^) 

This  agrees  with  Eq.  (5«3)' 
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The  total  ground  state  energy  in  zero-order  is  the  sum  of  the  sepa- 
rate nuclear  and  interelectron  energies.   Thus^  using  Eq.  (3.6)^ 

E^  =  -  Z^  -  (1  +  bjq^ 
or 

E^  =  -  Z^  +  Z  -  1  (6.7) 

in  viev  of  Eqs.  (6.3)  and  {6.6).      The  corresponding  zero-order  energy  using 
only  the  hydrogenie  functions  is  just 

E^  =  -  Z^  .  (6.8) 

o 

Table  I  presents  the  zero-order  energies  calculated  using  interelectron  func- 
tions, and  also  using  only  hydrogenie  functions.  These  results  are  given  in 
terms  of  the  percent  deviation, A,  from  the  true  value: 


^  =  100  ^^^  "  "^^  .  (6.9) 

expt 

It  is  instructive  to  compare  the  relative  magnitudes  of  Vp(0, r  )  and 
H'(0,r  ),  [see  (^.2)].   They  can  be  most  easily  pictured  by  placing  the  inner 
electron  at  the  origin.   Z~  V  (0,r  )  and  Z  "TI  (0,r  )  are  independent  of  Z  for 
large  Z  and  have  the  values 


(6.10) 


V2(0,r2) 
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-,   r       -r 
3e  -  e 

^   r,      -r 
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(6.11) 
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The  last  equation  can  be  obtained  from  the  confluent  hypergeometric  solution 
given  in  Appendix  I  or  by  solving  Eq.  (3-5)  directly^  with  b  =  -  Z^  for  large 
Z,  taking  into  account  the  singularity  at  the  origin.   Fig.  1  presents  a  com- 
parison of  V_  and  H'  as  given  by  the  above  equations.   The  small  magnitude  of 
the  perturbing  hami Itoni an,  H',  relative  to  the  Coulomb  potential,  ¥_,  is 
quite  evident. 


b.   Excited  states  of  helium 

An  excited  tvo-electron  atom  is  the  simplest  atomic  system  involving  more 
than  one  orbit.   Consider  one  electron  in  the  inner  I       shell  and  the  other 
electron  in  the  outer  m   shell.   The  inner  electron  shields  the  outer  elec- 
tron from  the  nuclear  charge  but  sees  the  full  nuclear  charge  itself.   The 
zero-order  wave  function  is 

♦0=5^(^1)4^(^2)  ^m(-12)-  (^-^2) 

Zr^   Zr^ 
At  large  distances  this  has  the  asymptotic  behavior,  exp(-  —r—  -   +  qr  ^). 

When  the  outer  electron  is  far  from  the  nucleus,  the  shielding  produced  by 

Zr^   (Z-Dr^ 
the  inner  electron  requires  the  behavior,  exp(-  —r—  -   ).   Therefore, 

q  =  l/m.  There  is  no  equation  to  express  shielding  of  the  inner  electron  be- 
cause there  are  no  electrons  in  still  lower  orbits  than  it.  The  perturbation 
H'  has  the  asymptotic  value  [see  Eq.  (6.5)  and  (5'5)J 


H-  ^  \  +  %   .  (6.15) 

2    m      m 


As  in  the  unexcited  case,  we  must  choose  b  =  -  Z. 
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The  total  zero-order  energy  nov  becomes 


^  I  m  '      m 


The  corresponding  zero-order  hydrogen! c  energy  is 


^  I  m 


The  deviation  of  these  energies  from  the  experimental  values  are  presented  in 
Table  I  for  the  single  excited  symmetric  states  of  helium. 

c.   Ground  state  of  neutral  atoms 

We  now  treat  the  ground  states  of  neutral  atoms  vithin  the  periodic  chart. 
The  atoms  are  treated  by  starting  with  the  innermost  shell  and  filling  the 
shells  successively.   Those  orbits  with  the  same  principal  quantum  number  are 
considered  equivalent.   Once  an  inner  shell  is  filled  no  further  change  occurs 
in  the  manner  of  treating  that  shell. 

The  first  shell  is  filled  with  helium.   The  contribution  E  ^  to  the  zero- 

ol 

order  energy  from  this  completed  shell  as  we  progress  along  the  periodic  chart 
changes  only  because  of  the  variation  in  Z.   This  contribution  has  already 
been  given  by  Eq.  (6.7)  as  -  Z  +  Z  -  1. 

Consider  the  electrons  in  the  second  principal  orbit  n  =  2.  As  one  of 
these  electrons  is  moved  to  infinity,  the  electrons  remaining  in  the  second 
and  first  orbit  shield  the  nucleus.   Since  the  q   are  all  identical,  Eq.  (5.2) 


yields  q^   =  1/2,   n  <  2. 
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The  perturbation  terms  linking  this  shell  to  lover  shells  must  vanish 
asymptotically.   The  requirement  is  fulfilled  by  Eq.  (5.5)-   There  are 
[p^CPo-  1)-  2j /2  interelectron  coordinates  linking  the  second  shell  to 
itself  or  lover  shells.   Therefore^  the  zero-order  contribution  to  the  ener- 
gy from  the  second  shell  is 


-■£'  fv"^  ,1  (•    ^2^  r^2'v^' 


\2-     '-^      \-^>    *%      [?-^)[-^^--ll-  (6.16) 


2 
When  this  shell  is  filled^  P  =  10  and  E   has  the  value  -Z  +  IIZ  -  55.   In  a 

similar  manner^  contributions  from  the  outer  shells  can  be  calculated  rapidly. 

For  instance,  the  third  shell  contributes  the  energy 


^o,=  #(^)^H-?)(^4^-^i 


vhlle  the   r       shell  has   the  energy 


\r  "     ~r      V  2""      y    "^  "T     V^  -  2^y  V        2 

r  r 


\/P   (P  -   1)        P     ,(P     .-   1)^ 
r  H     r^    r     _         r-1^    r-1 

(6.18) 


The  total  zero-order  energy  is  the  s\am  of  the  various  energies  of  the 
filled  or  partially  filled  shells.   The  corresponding  zero-order  contributions 
using  only  hydrogenic  functions  is,  of  course 
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2    fP  -  P    ) 
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The  energies  calciilated  with  the  above  formulas  are  given  in  Table  11^  in 
terms  of  their  deviation  from  the  experimental  values. 

7.   Discussion 

Table  I  presents  the  comparative  energy  deviations  calculated  in  zero- 
order  vlth  the  use  of  hydrogenic  functions  and  with  the  method  as  outlined 
in  this  paper  for  various  excited  and  ionized  two-electron  atoms.   As  noted 
in  Appendix  I,  our  method  is  not  applicable  to  the  case,  Z  =  1.   The  intro- 
duction of  interelectron  coordinates  together  with  the  proper  shielding 
requirements  produces  a  very  substantial  gain  in  accuracy.   This  is  illus- 
trated more  forcibly  in  Table  II  where  results  for  various  neutral  atoms 
are  given.   It  is  apparent  that  hydrogenic  functions  do  not  yield  satisfac- 
tory energies  for  neutral  atoms  in  zero-order. 

In  order  to  yield  improved  answers  with  hydrogenic  functions,  recourse 
must  be  made  to  first-order,  or  still  better,  to  variational  calculations 
with  a  consequent  large  increase  in  computational  effort  at  each  stage.   We 
also  present  in  Table  II,  these  first-order  and  variational  calculations  for 
the  hydrogenic  functions.   The  variational  results  are  available  and  were 
obtained  with  the  antisymmetrized  Morse  functions.'--'^  L  J   These  Morse  func- 
tions are  just  the  hydrogenic  wave  functions  in  which  the  constants,  such  as 
the  nuclear  charge,  are  treated  as  variable  parameters.   The  Morse  function 
calculations  are  representative  of  the  many  different  variational  calculations 
in  the  literature  aad,  furthermore,  have  been  extended  over  most  of  the  range 
for  which  energy  values  are  available.   The  first-order  energies  given  in  the 
table  were  obtained  by  using  the  formulas  given  in  reference  6 
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and  by  taking  the  wave  function  as  a  product  of  unsymmetrized  hydrogenic 
functions. 

Table  II  shows  that  the  zero-order  energies  obtained  with  the  intro- 
duction of  shielded  interelectron  functions  are  considerably  better  than 
the  first-order  hydrogenic  energies.   This  is  so  despite  the  fact  that  the 
zero-order  results  are  enormously  more  simple  to  obtain,  involving  only 
algebra  and  no  integration.   Indeed,  the  zero-order  energies  are  surpris- 
ingly close  to  the  variational  energies  in  view  of  the  number  of  variable 
parameters  used  in  these  latter  calculations. 

Another  widespread  and  relatively  simple  method  for  obtaining  the 

energies  of  neutral  atoms  is  the  Thomas-Fermi  calculation.   This  calculation 

is  based  on  a  variational  approach"-  -' ;    however,  because  of  the  mathematical 

approximations  introduced  within  the  theory^  -^ ,  accuracy  in  the  energy  values 

comparable  to  usual  variational  calculations  cannot  be  expected.   Table  II 

presents  the  deviations  calculated  by  the  usual  Thomas-Fermi  expression, 

E   =  O.769Z  '  .   These  results  compare  very  poorly  with  our  zero-order 
TF 

calculations  for  the  light  atoms  given  in  Table  II,  the  Thomas-Fermi  energy 
deviations  being  an  order  of  magnitude  larger  than  the  zero-order  deviations. 
The  experimental  energies  for  the  light  atoms  have  been  extended  somewhat  to 
heavier  atoms  by  Mayer  and  these  energies  have  been  reported  by  Cowan  and 
Ashkin'--'.   Fig.  2  compares  the  Thomas-Fermi  energies  and  our  zero-order 
energies  with  the  experimental  values  over  the  extended  range  of  Z  up  to 
those  values  of  Z  where  relativistic  corrections  become  important.   No  signi- 
ficant improvement  of  the  Thomas-Fermi  results  over  the  zero-order  calculations 
is  found  within  this  range. 
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Some  vords  should  be  said  here  about  the  applicability  of  the  type  of 
shielding  we  introduced  as  a  boundary  condition  upon  the  interelectron  wave 
functions.   This  discussion  may  help  us  to  understand  how  the  results  we 
obtained  can  be  adapted  to  other  types  of  quantum-mechanical  systems.   It 
should  be  remembered  that  the  outer  shells  were  neglected  in  treating  the 
asymptotic  shielding  of  any  orbit,  while  all  subshells  in  the  same  shell 
were  treated  identically.   Such  a  procedure  will  be  expected  to  work  well 
whenever  the  atomic  system  possesses  well  separated  shells  together  with 
an  approximately  central  field.   The  characteristic  of  the  Coulomb  poten- 
tial is  that  it  does,  indeed,  produce  well  separated  shells  and  the  central 
field  approximation  is  certainly  excellent  for  all  the  filled  shells  within 
the  atom.   Thus  our  approach  should  describe  very  well  the  behavior  of  the 
inner  shells  in  an  atom  and  these  shells,  of  course,  contribute  the  major 
part  to  the  total  energy  of  the  atom.   The  presence  of  p,d, ...  shells  out- 
side of  filled  shells  should  disturb  the  results  somewhat,  but  such  inaccu- 
racies will  disappear  once  the  shell  is  filled.   A  glance  at  Table  II  confirms 
these  expectations  in  a  general  way.   Note  the  extremely  good  results  once  a 
shell  is  completed  in  comparison  to  the  progressive  worsening  while  the  shell 
is  being  filled. 

A  final  note  will  be  made  concerning  the  applications  of  first-order  per- 
turbation theory  with. our  zero-order  wave  functions,  if   .      Because  of  the  inter- 
electron functions,  the  wave  functions,  \|/  ,  generated  by  allowing  the  individ- 
ual electrons  to  fill  all  of  the  possible  orbits  are  not  completely  orthogonal 
among  themselves.   This  is  so  because  H'  is  not  zero.   If  it  were,  we  would 
have  H  =  H  and  the  set,  ^     ,  would  be  orthogonal  simply  because  of  the  Her- 
mitian  properties  of  H.   However,  our  method  of  approach  has  been  aimed  at 
making  H'  small  compared  to  H  and  the  results  obtained  certainly  indicate 
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that  this   aim  has  been  accomplished.     We   can  then  be   fairly  confident  that 

the   set   of  fionctions,    \|f     ,    is   quasi-orthogonal  in  the   sense   that   (♦q-p^^qq^ 

is  much  smaller  than  either   (^     ,^      )    or   (itr       \|r      j.      This,   plus  the   fact 

^'op  op     ^  oq^  oq 

that  H  -^^  H,  will  allow  us  to  take  over  the  usual  first-order  perturbation 
o 

formulas  directly  for  use  with  the  set,  +  ,  with  the  expectation  that  any 
inaccuracies  will  affect  the  results  in  second-order  only. 
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Appendlx  1 


The  solution  to  the  equation  of  the  interelectron  wave  function 


2 

/-,-,,  cl  u    2  du    U    „     ^  ,  T  -,  , 

(1  +  b  — ^  +  -1 +  Eu  =  0  (I.l) 

^        ^  2   r  dr   r 
dr 


has  the  form 


"  -  ='^^.  (Vr§   '  T^  '  -  ^H  f^-^' 


where  q  =  ±  y/-E/(l+b)  and  F  is  the  confluent  hypergeometric  function.   At 
/large  values  of  r  and  for  real  q, e   ,   F  has  the  dependence  exp(+  |q|r)<. 
This  is  the  exponentially  increasing  solution  in  r  at  infinity  which  is  needed 
to  fulfill  the  shielding  requirement. 

For  the  hydrogen  negative  ion,  (b+lj  has  the  value  0.   The  solution  of 
Eq.  1  for  that  case  is  n(0)exp|_(Er  /4)-(r/2jJ.   This  does  not  have  the  expo- 
nential behavior  we  need  and  our  method  is  not  applicable  in  that  particular 
instance. 

It  is  worth  noting  that  the  integrals  necessary  to  obtain  the  atomic  ener- 
gies from  the  variational  formula  E  =  (\|f  .H>|r  j/(i^  ,^    )    can  be  carried  out  in 

o   o    o  o 


closed  form  for  two  electron  systems  when  (l.2)  is  used  as  the  interelectron 

[12] 

•S  '-       -*      VI 

[1]  ^  [2] 


ri2i 

wave  function.   For  these  systems  the  Hylleraas '-  -'  variables  are  introduced 


and  integrals  of  the  following  types  are  found' 


,J.,  =  /  e-(-'  /,   (-  ^^^   ,  .  ^  ,  -  ..)  .V  .   (X.5, 
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These   are   evaluated  by  the   formulas   given  by  Alder  and  Winther'-    -• 
The   final   result   is 


n-  TP    f     ^1  -1  -1  -2  -2  -1  -1     \       ,      ,   , 


[2(Z-ljJ  V  2(Z-lj        2(Z-1J        Z-1        Z-1        Z-1        Z-1. 

where  F^  is  the  generalized  hypergeometric  function  in  two  variables'-  -' . 

An  approximation  technique  was  given  in  I  for  evaluating  integrals 
containing  many  interelectron  functions.   For  greatest  accuracy,  the  tech- 
nique depended  upon  the  evaluation  of  integrals  of  the  type,  J  .   Eq.  (l.'+j 
thus  allows  this  approximation  technique  to  be  used  with  the  wave  functions 
of  this  paper. 
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Appendix  II 

Because  of  the  presence  of  the  interelectron  wave  functions  u. .  vhich 
vary  in  form  vith  changes  in  the  orbits  which  the  electrons  i  and  j  occupy, 
the  usual  Slater  determinant  cannot  be  constructed  in  order  to  represent  the 

ri5i 

antisymmetrization  of  \1«  .   However,  an  alternative  prescription'-  -J  suffices 
for  the  problem. 

Define  the  numbers  e  „     as  zero  when  any  two  subscripts  are  equal  and 
as  +1  or  -1  when  the  a,p...v  are  all  unequal  and  form  either  an  even  or  odd 
permutation  of  1,2, ...N.   For  our  purposes,  we  assume  that  the  ^  includes  both 
the  hydrogenic  functions  and  the  spin  functions »  We  label  these  ^'s  with  a 
Greek  subscript  and  since  the  form  of  the  u. .  depends  only  on  the  orbits  which 
they  connect  we  also  switch  the  subscripts  on  the  u  to  the  Greek  subscripts 
corresponding  to  the  connected  orbits.   Let  1,2,...   represent  the  nuclear  and 
spin  coordinates  of  the  first,  second,  etc.  electron  and  12,13,...   the  cor- 
responding interelectron  coordinates. 

The  antisymmetrized  zero-order  wave  function  is 

*c  =  ^a^...v  t^^^  lp(2)...  ^N)  u^(l2)  u^^(l5J...  u^^(MN)  •    (ll.l) 


Any  subscript  which  appears  more  than  once   indicates   a  summation  over  all  values 
of  that  subscript. 

That   \|r     is   truly  antisymmetric  requires  the   reasonable   assumption  that 


u^    (FD)    =  u^    (DF)    =  u  ^(FD)    =  u  -(DF)   .  (ll-2j 


This  is  in  conformity  with  the  ideas  set  forth  in  this  paper. 
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Table  I 


SYS'i'KM 

E    (a.u.j 
exp  ^ 

0 

A 
0 

Helium-like 

He 

-  2.90)+ 

-37.9 

-  3.3 

Li-^ 

-  7.29 

-23.5 

+  i+.o 

Be^-^ 

-15.67 

-17.0 

+  k.9 

b5-^ 

-22.0^ 

-13. i+ 

+  h.7 

C^^ 

-32.4 

-11.1 

+  i+.3 

IT  5-^ 

-1+4.8 

-  9.h 

+  4.0 

0^^ 

-59.2 

-   8.1 

+  3.7 

F^^ 

-75.5 

-  7.2 

.3.^. 

Excited  Helium 

o 

-  2.162 

-15.6 

-  4.1 

-  2.065 

-  7.5 

-  2.2 

o 

-  2.059 

-  ^.3 

-  1.2 

Table  I.   Comparison  of  percentage  deviations,  A  ,  A  of  calculated 
zero-order  energies  from  the  experimental  energies  using  pure  hydrogenic 
functions  and  using  the  method  of  this  paper,  respectively.   The  experi- 
mental values  for  excited  heliinn  are  an  average  of  the  singlet  and  tri- 
plet values. 


Table   11 


SYSTEM 

E   (a.u.) 

0 

4 

V 

^TF 

c> 

He  ■''S 

0 

-   2.90^+ 

-57.9 

+  5.5 

+l.9(lj 

-55.5 

~     1  •      1 

Li  ^S 

0 

-  7.^9 

-55.2 

+  6.7 

+0.9(5) 

-55.5 

+1.5 

Be  -^S 
o 

-  1^1.68 

-56.2 

+  7.8 

+0.8(4) 

-53.0 

+1.2 

B  2p 
o 

-  2*^.67 

-39.5 

+  9.8 

+0.6(4; 

-33.5 

-0.5 

C  ^P 
o 

-  57.9 

-42.6 

+11.5 

+0.6(4) 

-52.7 

-1.6 

N   S^ 

-  5^.6 

-45.7 

+12.9 

+0.6(4) 

-52.0 

-2.7 

0  ^Pg 

-  75.2 

-49.0 

+14.2 

+0.8(4) 

-30.9 

-5.7 

2 
F   P^ 

-  99.9 

-52.0 

+15.2 

+1.0(4) 

-29.7 

-4.6 

Mg\ 

-200.3 

-52.1 

- 

- 

-26.6 

-0.5 

2 

Al  P 

0 

-242.7 

-51.0 

- 

- 

-25.9 

+0.2 

TT  TT  TT 

Table   II.      Percentage   deviations    calculated  for  neutral   atoms.      A     ,    A,    ,    A 
are   the   zero-order_,    first-order  and  variational  results   using  hydrogenic   func- 

tions.      The  n\imbers   in  parenthesis   after  A     are  the  nimber  of  variable  para- 

^       TF 
meters   used  in  the  variational  calculation.      A       and  A     are  the  Thomas-Fermi 

o 

results  and  the  zero-order  results  on  this  paper,  respectively.   The  experi- 
mental energies  are  obtained  from  C.E.Moore  L  J  g^  extended  somewhat  by 
Cowan  and  Ashkin  "-  -'  . 
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